Thermal conductivity in B- and C- phase of UPt 3 
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Although the superconductivity in UPt3 is one of the most 
well studied, there are still lingering questions about the nodal 
directions in the B and C phase in the presence of a mag- 
netic field. Limiting ourselves to the low temperature regime 
(T^C A(0)), we study the magnetothermal conductivity with 
in semiclassical approximation using Volovik's approach. The 
angular dependence of the magnetothermal conductivity for 
an arbitrary field direction is calculated for both phases. We 
show that cusps in the polar angle dependence appear in B 
and C phases which are due to the polar point nodes. 

PACS numbers: 74.20.Rp, 74.25.Fy, 74.70.Tx 



I. INTRODUCTION 

Perhaps due to the exceptional appearance of two su- 
perconducting(sc) transitions with two associated criti- 
cal field curves, UPt3 is one of the most well studied 
Heavy Fermion superconductors. The history of this sub- 
ject was described in Ref.El. The T linear dependence of 
the low temperature thermal conductivity both paral- 
lel to the c axis and the a axis was in famCirQf the Fi2u 
rather than Ei g - type sc order parameter&LKi Further 
Pt- NMR experiments confirmed the triplet nature of the 
superconductivity and the details of the triplet gap fuai 
tion d (k) n A, B, and C- phases were determinecE'l 
This assignment of the spin configuration is consistent 
with the weak spin-orbit coupling limittfflj also assumed 
here. H C 2 Along the c axis requires a strong Pauli lim- 
iting effectEI. This favors strong spin orbit coupling but 
then H C 2 results are inconsistent with the NMR data. We 
assume that d (k)=A(k)z to be weakly pinned along c. 
Then the E2« gap functions of for the low temperature 
B and C phases of UPt3 are given by 

q 

B : A(0,p) = -v / 3Acosi?sin 2 ??exp(±2^) 

C : A(<8,<p) = ^v / 3Acos t ?sin 2 i?cos(2v3) (1) 

Here ip are the polar and azimuthal angles of k . The 
angular dependence of A($, ip) is shown in Fig. |l|. In the 
B- phase the poles are second order node points and the 
equator is a node line. In the C- phase two additional 
vertical nodal lines appear at 45 degrees away from the 
vertical plane containing H . The thermal conductivity in 
the vortex phase for field along the symmetry axis also 



decided in favor of the E2,j-state and against the Ei 3 - 
statcEJ, this is reviewed incil. The angular—dependent 
thermal|-cpnductivity has been measured ir£3 and ana- 
lyzed inIL 3 ! Unfortunately the experiment was done for 
T>0.3 K which is not sufficiently low to determine node 
structures. Recently we have been stu dying the , th ermal 
conductivity in nodal superconductorsata-MO'lla'tj. In 
particular we shall restrict to the low temperature limit 
T< vVeH < A(0) and the superclean limit (TA)i < 
wV eH, where v = (v a v c )^ , T is the scattering rate and 
v a ,c are the anisotropic Fermi velocities. The condition 
v\J eH <C A(0) can be satisfied in the B- phase while 
vV eH < A(0) in the C- phase since the latter appears 
only for H> 0.6 T and H> 1.2 T for field along a and 
c respectively. Our results in the C- phase may there- 
fore only have qualitative significance. Restriction to the 
superclean limit does not influence the main conclusions 
on the connection between node topology and magne- 
tothermal conductivity. We first examine the quasipar- 
ticle DOS and the thermodynamic properties of the B 
and C phase of UPt3 in the presence of a magnetic field 
with arbitrary orientation at low temperatures. Then we 
study the thermal conductivity in the B and C phase 
in the low temperature regime which provides clear ev- 
idence for the nodal positions in A(k) . In this regime 
the influence of the AF order in UPt3 may be neglected 
contrary to the A- phase regime. 



II. DENSITY OF STATES AND 
THERMODYNAMICS 

The quasiparticle DOS in the vortex state of UPt 3 is 
given by 



g(Q) = Re( 



C - ix 



V^Cc - ix) 2 + f 2 



(2) 



Where the brackets denote both Fermi surface and vor- 
tex lattice average. With the form factor fs(z)= 3 -^(l — 
z 2 )z and f c (z,z')=(2z' 2 -l)f B (z) (z=cosi9, z'=cos<p) for B 
and C- phase respectively the average may be computed 
and one obtains 
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where x— ^ v A q I is the normalized Doppler shift of 
quasiparticle energies, Co = lim^^o Im(Q/A) and u> is 
the renormalized quasiparticle energy in the presence of 
impurity scatteringESEJ. Furthermore v is the quasipar- 
ticle velocity and 2q is the, pair momentum around the 
vortex. Following VolovikEJ wc obtain 



(x) B = -v^I B (9) 

7T A 



2 2 
(.xln(-)) c = -v— r- 

X 7T A 



Ib(8) = ol sin 
Icifi) =I B (0)+I(6) 
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7T 
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Here E(sin#) is the complete elliptic integral of the 
second kind. Furthermore we used the definition 



1 
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+ cos( 



— or sm 

tan 

v 2 



2a sin 9 



L (atan<9))} 



(5) 



/(*) = 



cos 1 (a sin 6*) for yfasm(8) < 1 
cosh - (a sin 9) for i/asm(9) > 1 



where a = v c /v a is the anisotropy of Fermi velocities 
and 9 is the polar angle of H with respect to c axis. The 
averages in Eq. ^ are evaluated by replacing the k space 
integration by a summation over nodal positimi^and then 
integrating out the superfluid velocity fieldtjo. We as- 
sume a square vortex lattice for simplicity a hexagonal 
lattice would lead to an additional numerical factor 0.93. 

Substituting Eq. | into Eq. | we obtain 



9b 



9c 
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7T 
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A 



-Ib(9) 
W Ic(0)H^ 



(6) 




FIG. 1. Spherical plots of |A(tf, ip)\ for B phase (left) and C 
phase (right) . For the C- phase two additional vertical nodal 
planes appear an angle j away from the vertical plane which 
contains H assumed at d> = 5- 



In addition we have Co=£[<7(0)] 1 . The low temper- 
ature specific heat, the spin susceptibility etc. are given 
by 



Cs Xs 1 Ps(H) 
InT xn Ps{0) 



(7) 



where p s is the superfluid density. The 9- dependence 
of g(0) for B- and C- phases is shown in Figs. || and || 
respectively. 



III. THERMAL CONDUCTIVITY IN THE 
VORTEX PHASE 

The thermal conductivity tensor in the vortex phase 
depends on the angles (9, 0) of H due to the angle de- 
pendence of the Doppler shift energy Ax. Here cj> is the 
azimuthal angle between H and the direction of the heat 
current jq in the ab- plane. Following)^ we obtain for 
the B- phase 



Kzz 
Kri 
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F X B X { 



l V -^(eH)I B (9)F^(9) 
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FIG. 2. Polar field angle dependence of Is(6 l ) and 
Ib(0)F % £(9, 45) (ij= xx, zz, xy) which determine the 9- de- 
pendence of DOS g(0), thermal conductivities (k xx , k zz )) and 
thermal Hall coefficient (n xy ) respectively. 

The thermal Hall coefficient in the B- phase is obtained 



2 



F x B y { 



3A 2 
2 sin(20) 
7T 3sin 2 

[(2 - sin 2 9)E(sm 9) - 2 cos 2 01f (sin (9)] 



(9) 



The 9- and 0- angle dependences of (ij=xx,zz,xy) 
in the B- phase are shown in Fig. || and Fig. For heat 
current along c (k zz ) no </>- dependence appears. 

In the limit 9 — ^, Is( : |)=a + — and then 



Kxi ~ icos(20)); 

7T O 



— sin(20) (10) 



The maximum in k xx (90, 4>) at cf>= ±90 occurs for heat 
current _L H when the Doppler shift gives rise to the 
largest quasiparticle DOS parallel to the heat current and 
we have k xx {4> = 90) / k xx (4> = 0)=2. 




FIG. 3. Polar field angle dependence of lc(9) and 
I O (0)F«(0,45) (ij=xx, zz, xy) 

Now we consider the C- phase. Again k zz does not 
exhibit a <fi- dependence. The C-phase according to Eq. |l| 
has two additional perpendicular node lines. Rotating 
the field at a given 9 around c (i.e. changing <f) will 
lead to a co- rotation of these node lines such that the 
vertical plane containing H always stays at half angle 
between the.-two perpendicular planes of the node lines 
parallel to cEI(see inset of Fig. ||). Consequently k zz (9) 
will again be independent of <fi while k xx {9, (/>) depends 
on both field angles. We find for heat current along c: 



K>n D A VVeH 



F^ z {9) =asin0 + - / dz\z\ [-(1 + cos 2 0)(1 - z 2 ) + 
7T J-i 2 

(a 2 sin 2 (9z 2 + V2 sin (9 cos 6»z(l - z 2 )?] * (11) 
On the other hand we obtain for heat current along a: 



k xx 



V2(l 



= F% X { 



XX / 

C \ 



COS 



'ln 2 (^^) (12) 



As is readily seen k zz depends only on 9, while n xx 
depends both on 9 and </>. Both angular dependences are 
shown in Fig. || and Fig. |5| respectively. 




FIG. 4. Azimuthal field angle dependence of l B (9)F% X (6, <j>) 
for various 9. For 9 — ►() the (j> dependence is suppressed com- 
pletely. The zz- component is always <j>- independent. 

In both B and C phase «y (ij= xx,zz) and also the 
specific heat which is determined by 1b,c exhibit clear 
cusps at the poles 9=0 (and 9 = n) caused by the con- 
tributions from the respective second order node points 
which are present in both B and C phase. This is very 
similar to what has been recently observed and analyzed 
in YNi2B2Cc3. There the second order node points lie 
along the equator and hence the cusps appear as func- 
tion of <j>. The most significant difference in the B- and 
C- phase results can be seen in the behaviour of k zz for 
small polar angle 9. While in the B- phase it approaches 
zero it remains finite in the C- phase. Furthermore the 
anisotropy for 9 = ^ is considerably smaller 
in the C- phase as compared to the B-phase. In both B- 
and C- phase the xx- component exhibts nonmonotonic 
behaviour as function of 9. 

The thermal Hall coefficient in the C- phase reads 



(eH)I c (9)F* 



(13) 



where Y X J( 



) = Fg V (9,(f)) holds because the con- 
tributions from perpendicular node lines with H lying at 
half angle in between cancel and so as in the B- phase one 
is left with polar and equatorial contributions to the ther- 
mal Hall constant. For numerical calculations we used 
the anisotropy ratio a = ^f- =1.643. It can be directly 
obtained from experimental anisotropies of thermal and 
electrical conductivities^ — = -2^=2.7 which are equal to 

~,2 



3 




FIG. 5. Azimuthal field angle dependence of l c (d)Fc x (d, <f>) 
for various 8. The zz- component is again 4>- independent. 
The inset shows the field and node geometry in the C- phase. 
The node lines along c are lying in the planes perpendicular 
to ab (dashed lines) which are mutually orthogonal. The field 
H lies at half angle ( -f ) in between forming an angle <f> with the 
heat current jq along the a- axis. This geometry is preserved 
for any 4> because the nodal planes co- rotate with the field. 



IV. CONCLUDING REMARKS 

We have foun d that in UPt3 at low temperatures 
(i.e. T<C vVeH) and in the superclean limit ((rA)5 <g; 
v\feH), the thermal conductivity exhibits clear angular 
dependence which will help to identify the nodal direc- 
tions in A(k) and to verify the predictions of the com- 
monly discussed E2 U model for the gap function in B 
and C phases. Most significantly we predict that i) cusps 
appear in the thermal conductivity and specific heat for 
9=0 (and 9 = ir) due to the polar point nodes of UPt3 in 
both B- and C- phase, ii) in the C- phase for heat cur- 
rent along the nodal direction (<f> = ?) a finite k zz (9, j) 
occurs even for 9 — >0 which is caused by the contribution 
from the additional perpendicular node regions. For the 
B- phase this contribution vanishes, iii) the k xx — k zz 
anistropy for H in the ab- plane (6 = ?) is considerably 
larger in the B- phase as compared to C- phase. This 
is again caused by the perpendicular node lines which 
contribute and enhance k xx only for the C- phase. 

We hope that these different features in the field- angle 
dependent thermal conductivity tensor above and below 
the critical field of the B-C transition will resolve a part 
of the remaining controversy surrounding A(k) in UPt3 . 
We recall that thermal conductivity experiments have 
been very useful to identify A(k) in unconventional su- 
perconductors. For example, Izawa et ai-have succedcd ta 
identify A(k) in S^RuOja CeCoIn 5 a, organic saltsEa 
and more recently YN12B2CO. Indeed the thermal con- 
ductivity appears to provide a unique window to access 



the nodal structure in unconventional superconductors. 
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